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ABSTRACT 

Far-ultraviolet (900- 1200 A) spectral synthesis of nine giant extragalactic H II regions in M33 and 
M 101 is performed to study their massive stellar content. Several parameters are quantified, predicted, 
and compared to the literature: age, stellar mass, IMF slope, number of O-type and Wolf-Rayet stars, 
Ha and 5500A continuum fluxes. The results of this particular technique are consistent with other 
methods and observations. This work shows that a total stellar mass of a few 10 3 M Q is needed 
to populate the IMF bins well enough at high masses to obtain accurate results from the spectral 
synthesis technique in the far- ultraviolet. A flat IMF slope seems to characterize better the stellar 
line profiles of these objects, which is likely the first sign of a small number statistic effect on the IMF. 
Finally, the H II region NGC 5461 is identified as a good candidate for hosting a second generation of 
stars, not yet seen at far-ultraviolet wavelengths. 

Subject headings: stars: early-type — galaxies: stellar content — galaxies: individual (M33, M101) 
- ultraviolet: stars 



1. INTRODUCTION 

Giant extragalactic H II regions (GEHR) are impor- 
tant sites of star formation. They are small scale ex- 
amples of extreme sites of star formation such as lo- 
cal and distant starburst galaxies. Like starburst re- 
gions, they contain several distinct st ar clusters (e.g. 
iMeurer et al] ri995: H unter et all Il996j) that can inter- 
act with each other to potentially enhance or slow 
down the star formation processes (see review by iTanl 
I2005L and references therein). They produce the most 
massive stellar types known (O, B, and Wolf-Rayet) 
that have the potential to transform the morpholog- 
ical and che mical aspects of galaxies through their 
feedback (e.e. iHeckman et all 119901 iMartin et alJl2002l 
iTremonti et all 120041 iCalzetti et all 12004 among oth- 
ers). Most GEHR are rec ent and quasi- instant a neous 
events of star- format ion jMa^JHe^s^^FCunib 119911 
119991 iSchaerer et all 119991: iStasihska fc Schaererl Il999|) . 
as seems to be the case for a starburst, in the sense that 
most of t heir massive st ars seems to form within less than 
2-3 Myr l)Pellerinll200^) . 

Evolutionary synthesis is a powerful tool to 
study stellar popula ti ons in various enviro n ments 
fe.g. IWorthev et all 119941 ILeitherer et all 119991: 
IBruzual fc Charlotl 120031 IRobert et all 12003ft . The 
main goal of evolutionary synthesis is to deduce the 
global properties of spatially unresolved stellar popula- 
tions such as their averaged age, mass, and metallicity. 
The development of evolutionary synthesis codes in 
the past decade has con siderably improved our knowl- 
edge of galaxies (e.g. iGonzalez Delgado et all 119991 
ILeitherer et alJl2001t IChandar et alJl2003l among many 
others). With the recent (and coming) generation of 
large telescopes such as Keck, Gemini, JWST, and 
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ALMA, this technique will be very useful for our under- 
standing of very distant galaxies and of their evolution 
through cosmic time. 

Nearby GEHR are excellent candidates to test the ac- 
curacy of the evolutionary synthesis technique. GEHR 
like those found in M 33 are close enough to resolve indi- 
vidual stars and to compare their detailed stellar content 
with what is deduced from synthesis of integrated spec- 
tra. In this work, a detailed study of the massive stellar 
content of several GEHR observed in the far-ultraviolet 
(900-1200A; FUV) is presented. The study is based on 
the spectral synthesis code LavalSB and i ts recent em- 
pirica l spectral library in the FUV range ({Robert et alJ 
2003). The synthesis of GEHR observed in M33 and 
M 101 will be compared, when possible, to previous works 
detailing their resolved stellar content. 

The following section presents a summary of the data 
processing. Section |3 describes the evolutionary synthe- 
sis code LavalSB used in this work. The synthesis results 
for each GEHR are detailed in 21 and compared with 
previous works at various wavelengths. A discussion of 
specific results is presented in section 5 and the main 
results are summarized in §6. 

2. FUSE DATA AND REDUCTION 

FUV spectrograms of nine GEHR were obtained by the 
Far Ultraviolet Sp ectroscopic Explorer (FUSE) telescope 
l)Moos et al.ll200f)ft for various projects. Most data were 
obtained through the largest aperture (LWRS; 30"x30") 
while some spectrograms were obtained using a smaller 
aperture (MDRS; 4"x20"). Aperture locations are dis- 
played in Figure 1. A general description of the FUSE 
data is reported in Table 1. Data were gathered from 
the MAST 1 public archives. The data were processed 
with the calfuse pipeline v2.4.2. This version cor- 
rects for Doppler shift induced by the heliocentric mo- 
tion of Earth, event bursts, the walk problem, grating 
thermal shifts, bad pixels, background noise, distortions, 

1 M ultimission Archiv e at Space Telescope Science Institute; 
|http: / / archive. stsci.eduTl - 
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and astigmatism. More information relative to calf use 
is available electronically 2 . The output from calf use 
comprises eight segment spectrograms for each exposure 
that correspond to the eight optical paths of the instru- 
ment. Each segment covers a different wav elength range, 
with some of them overlapping (see fig. 2 of S ahnow et al.l 
2000). First, for each segment, each exposure was com- 
bined with a statistical weight based on exposure time. 
Then the segments that cover the same wavelength re- 
gions (roughly 900-1000A, 1000-llOOA, and 1100-1200A) 
were averaged with weights based on their signal-to-noise 
ratios. Finally, the spectrograms of each wavelength 
range were simply coadded to obtain one spectrogram 
covering the entire 905-1187A range. The spectrograms 
were then smoothed by a factor of 20 using the IRAF 3 
boxcar task, corresponding to aresolution of about 0.13A. 
This last step increases the signal-to-noise ratio without 
affecting the stellar line profiles. The spectrograms were 
corrected for rcdshift. Reddening correction will be dis- 
cussed in section 0J 

3. STELLAR POPULATION MODELING IN THE FUV 

A first wor k of spectral synthesis below 1200A has been 
made by iGonzalez Deleado. Leitherer. fc Heckmanl 
(|1997t) for the O Vi+Ly/3+C II feature. The stellar 
library was based on Copernicus and the Hopkins 
Ultraviolet Telescope (HUT) with a spectral resolution 
of 0.2A. Their work clearly showed that the line profile 
was sensitive to the age of a stellar population. A 
new FUV spectral library, based on FUSE data, has 
recently been added to the spectral synthesis code 
LavalSB ilh'obert et all 120031) . This code has been 
proven to be ve ry powerful for young stellar populations 
(jPellerinl [2004) and will be used in the present work 
to deduce the global properties of massive stars in 
GEHR from their integrated FU V light. LavalSB is a 
parallel version of Starburst99 (Leithe rer et al.l ll999). 
It uses the evolutionary tracks of the Geneva group 
(Schaller et al. 1992; Schaerer et al. 1993a, 1993b; 
Charbonnel et al. 1993; Meynet et al. 1994). The 
stellar population follows a mass distribution based on a 
chosen stellar initial mass function (IMF) and mode of 
star formation (instantaneous or continuous). Individual 
stellar parameters are used to assign the corresponding 
normalized empirical spec trogram from the FUV library 
based on relations from iSchmidt-Kalerl l)1982H . The 
normalized library spectrograms are flux calibr ated 
using stellar atmos phere models of iKuruczl l)1992|) for 
normal stars, and of lSchmu tz_et_al.l (1992) for stars with 
extended envelopes. The IKuruca 1 1992 ) spectra have 
been fitted using a Legendre function to remove their 
low resolution spectral features in order to avoid any 
confusion with empirical stellar lines from the spectral 
library. The FUSE stellar library covers from 1003.1 
to 1182. 678A with a dispersion of 0.127A. The library 
metallicities corresponds to the evolutionary tracks of 
LavalSB, e.g. Z<p fo r Galac tic stars (12+log[0/H]=8.7; 
lAllende Prieto et alJ 120011). 0.4 for LMC stars 
d2+logfO/H1=8.3: IRussell fc DopitallT99l . and 0.1 Z Q 

2 http:/ /fuse. pha.jhu.edu/analysis/calfuse. html 

3 Image Reduction and Analysis Facility, supported by NOAO 
and operate d by AURA u nder cooperative agreement with the 
NSF; |http: //iraf.noao.edu/| . 



for S MC stars (12+log[O/H]=8.0; IRussell & Dopltal 
1992). 

The most useful stellar indicators in the FUV are the 
C III blend multiplet centered at 1175.6A, and the P V 
doublet at 1118.0 and 1128.0A. The profiles of these lines 
show strong variations with age and metallicity of the 
population, depending on what spectral types dominate 
in flux. Significant, but more subtle, changes also appear 
with different IMF parameters. At shorter wavelengths, 
the O vi AA1031.9, 1037.6 and the S iv AA1062.7, 1073.0, 
1073.5 line profiles show variations with age and metal- 
licity, and possibly with IMF. However, the empirical 
stellar library used in LavalSB contain stars for which 
these diagnostic lines are contaminated by interstellar 
features from Galactic H2 and other atomic transitions. 
Consequently, stellar lines of O vi and S iv will not be 
used in the present work since C III and P V lines alone 
will provide more accurate results. An extensive identi- 
fication of stellar and interstell ar lines contained within 
the FUSE range can be found in lRobert et~aTl (|2003t) and 
iPellerin et all 1)20021) . 

To establish the characteristics of an integrated stellar 
population in the FUV, the FUSE spectrogram is first 
normalized and the stellar indicators of C III Al 175.6 and 
P V AA1118.0, 1128.0 are compared to the models. The 
best-fit model is chosen both by eye and by performing 
a x 2 fit- This first step provides information on the age, 
the metallicity, and the IMF parameters of the popula- 
tion. A standard IMF is defined here as having a slope 
a=2.35 a mass range from 1 to 100 M Q . Once the age 
and metallicity of the stellar population are estimated 
from the normalized FUV spectrogram, the extinction 
is then evaluated by comparing the observed continuum 
slope of the flux calibrated data t o the one of the best-fit 
model. The theoretical law from lWitt k Gordonl (|2000f) 
for a clumpy dust shell distribution with an optical depth 
of 1.5 in the V band is used to derive the internal extinc- 
tion E(B- V)j. The Galac tic extinction is corrected using 
the law of iSeatoiil l)1979|) . Finally, the stellar mass in- 
volved in the system is estimated from the unreddened 
flux level. 

Uncertainties related to the line profile fitting are de- 
termined by comparing the different sets of models at 
a given metallicity. Since LavalSB covers only specific 
values (0.1 Z Q , 0.4 Z , Z Q , and 2Z Q ), it is not possible 
at this point to evaluate the full age range that could fit 
the data. The jumps in metallicity are quite large so the 
synthetic spectra from the next metallicity value do not 
always reproduce the observed line profiles and cannot 
give clues on the age range. Consequently, the age un- 
certainties given in the present work are underestimated 
and do not take into account the possibility that the data 
can be fitted using a slightly different age at a slightly 
different metallicity. 

For a given model, the primary source of error in the es- 
timation of stellar masses and predicted fluxes is usually 
the FUV flux uncertainty from FUSE, which is usually 
around 10%. However, in some cases, the age uncertainty 
gives a larger error bar than the FUSE uncertainty. In 
every case, the largest uncertainty is given. Also, the 
IMF slope used to calculate the total stellar mass affects 
the uncertainty on masses and predicted fluxes. How- 
ever, these uncertainties are not explicitly included for 
the best fit model uncertainties. Where possible, param- 
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eters of other good-fit models are given to better evaluate 
the full uncertainties. 

4. MASSIVE STELLAR CONTENT OF GEHR: THE FUV 
POINT OF VIEW 

4.1. NGC604 

NGC604 is a well-known GEHR within the Local 
Group galaxy M33. Several studies found and con- 
firmed the presence of very massive O, B, and Wolf- 
Rave t (WR) stars llVilchez et, all 119881 iDrissen et alJ 
119931 | H unter et al.l 119 961 IGonzalez Delgado fe Perezl 
boodlBruhweiler et all2003HMai'z-Apellaniz et aU 2004). 
At least fou r distinct starcluste rs have been identified in 
this object l|Hunter et al.l ll9961. 

The FUSE spectrogram of NGC 604 obtained through 
the LWRS aperture is shown in Figure 2a. This aperture 
corresponds to a physical size of 123x123 pc 2 (l"=4.1p c 
at 840 kpc; see also Fig.l of iLebouteiller et all [2 005). 
The a perture includes the Cluster A from iHunter et all 
(1996), but not the entire H II region. The spectrogram 
has a very good signal-to-noise ratio (S/N) of 20 between 
1155 and 1165A that allows to perform a good synthe- 
sis with details on the IMF slope. The C III line profile 
shows a large absorption feature in its blue wing, indi- 
cating the presence of evolved late-type O stars. The 
P V doublet also displays PCygni line profiles typical 
of massive stars with strong winds. The line depths of 
C III and P V suggest a sub-solar metallicity for the 
stars. Continuous burst models have to be excluded 
since they produce stellar lines with too faint P Cygni 
profiles. To obtain a good fit, especially for the C III 
line profile, a flatter IMF with a slope a=1.5-2.2 is bet- 
ter, while a standard IMF with a=2.35 could also fit. 
The best fits are obtained for models having a=1.5 and 
an age of 3.9±0.1Myr for 0.1 Z Q and 3.3±0.1Myr for 
a 0.4 Z Q metallicity. If a(IMF)=2.35, then the best-fit 
ages are a little lower with 3.5±0.3Myr for 0.1 Z© models 
and 3.0 Myr at 0.4 Z©. The solution is not unique since 
there is a degeneracy in the line depth for the models at 
sub-solar metallicities when the P Cygni profiles are well 
d eveloped. 

IGonzalez Deleado fe Perezl ((2000) performed a de- 
tailed study of NGC 604 using IUE spectrograms 
(9.5" x 22" aperture), optical ground-based data, and Ha 
images from the HST, to fully describe this GEHR. From 
the H Balmer and He I absorption lines, they deduced 
an age between 3 and 4 Myr for the stellar population 
with a standard IMF or flatter. A continuous burst 
cannot fit their emission line ratios. Their IUE spec- 
trograms revealed a population of 3-5 Myr (better fit at 
3 Myr ) with an IMF slope flatter than 3.3. iVilchez et alJ 
( 1988) studied in detail the chemical abundances in M 33 
from nebular lines. They measured an oxygen abundance 
12+log[0/H]=8.51 for NGC 604. All these results are 
fully consistent with the FUV line profile synthesis. The 
best-fit model parameters are reported in Table 2, to- 
gether with the other good-fit models. Note that here- 
after, calculations using the models at 0.4 Z© are fa vored 
based on the metallicity from IVilchez et alTl|1988fl . 

Adopting an instantaneous burst model of 3.3 Myr at 
0.4 Z© with an IMF slope a=1.5, the observed FUV con- 
tinuum slope suggests no significant internal extinction 
E(B-V)j. No internal extinction is needed if a Galactic 
correction of 0.02 is applied, and E(B-V)i=0.03 is cal- 



culated if no Galactic extinction is applied. Using an 
IMF truncated between 1 and 100 M©, the FUV flux 
level leads to a stellar mass of (7±2) x 10 3 M© within 
the LWRS aperture. Using an IMF slope of 2.35, 
the calculated stellar mass i s rath er (1.4±0.3) x 10 4 M©. 
IGonzalez Delgado fc Perezl l|2000|) obtained a E(B-V); 
of 0.1 based on their IUE spectrograms. They also es- 
timate a stellar mass of 0.1-2xl0 5 M©. IHunter et all 
( 1996) found, based on optical HST images, an extinction 
value of 0.08 for Cluster A contained within the LWRS 
aperture. Their extinction and mass values are slightly 
higher than those from the FUSE data. 

From the stellar population described above with 
a=1.5, several physical parameters can be deduced 
and compared (see Table 3). First, such a popula- 
tion would theoretically lead to an unreddened Ha flux 
of (2±1) x 10~ n ergs s _1 cm~ 2 A -1 , and a continuum 
level at 5500A around (3±1) x 10~ 12 ergs s _1 cm~ 2 A^ 1 . 
Changing the IMF does not change these numbers sig- 
nificantly. Ha fluxes of 4.0 and 3.3 xl0~ n ergs s" 1 
cm~ 2 A -1 have been measured from HST and gro und- 
based image s bv IGonzalez Delgado & Perezl l)2000D and 
iBosch et al.l l)2002lL respectively. Those values are 
slightly above the FUV predicted values. The differences 
in Ha fluxes are consistent with the differences in stellar 
masses. 

Ac cording to the Ha+UV images fr om HST (see Fig. 
2 of IGonzalez Delgado fe Perezl 12000^ . several massive 
stars from Cluster A are co-spatial with the nebular emis- 
sion. These stars are good candidates to higher extinc- 
tion and it is likely that their contribution to the FUV 
flux is significantly lower than at longer wavelengths 
(even at ~1500A) and then partly explain the differ- 
ences observed in extinction v alues at various wave length 
ranges. Furthermore, Fig. 2 o flHunter et al.l 1^)96) shows 
that Cluster B and C contribute significantly to the 
nebular emission of NGC 604, but they are not taken 
into account in the total stellar mass derived from FUV 
since they are not included within the FUSE aperture . 
Also, a detailed study from iMaiz- Apellaniz et aLl (2004) 
revealed an extremely complex gas/dust geometry for 
which around 27% of the ionizing photons might be miss- 
ing in NGC 604 due to attenuation. In addition to the 
aperture effect, this obviously contributes to create a dis- 
crepancy between the predicted and observed values in 
the stellar mass and other fluxes parameters. 

The FUV synthesis of a 3.3 Myr population with a=1.5 
at 0.4 Z© predicts that about 9 WR stars (3 WN and 6 
WC) should be present in NGC 604. IDrissen et all ltl99l 
obtained ground-based and HST-WF/PC1 images and 
identified 12 WR or Of candidates, slightly more than 
the LavalSB predictions. More recently Drissen et al. 
(2005, in preparation) confirms that there are at least 6 
WN and 2 WC stars among them. This WC-to-WN num- 
ber ratio is not consistent with LavalSB (or Starburst99 
neither). To obtain WC/WN~l/3, both models propose 
an age around 4.5-4.7 Myr for the population. However, 
LavalSB do not include the effect of rotation in evolu- 
tionary tracks. By including rotation in the models, the 
result will be to extend the duration of the WR phase 
and to increase considerably the number of WN stars, 
which will fit better the observations (G. A. Vazquez 
2005, private communication). 
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Also, for the population synthesized above for 
NGC 604, LavalSB predicts that 90±?g O- type stars (of 
all spe ctral types still present at this age). iHunter et alJ 
(1996) estimate from HST/WFPC2 images that about 
190 stars brighter than 09.5 V are present in NGC 604, 
which is h igher than the FUV estimation. However, the 
number of IHunter et al.1 l)1996|) may include some B su- 
pergiants. If we use an IMF slope of 2.35, the model then 
predicts roughly the same number of O-type stars but no 
WR stars (or very few) at 3.0 My r, which is in dis agree- 
ment with the observations of lDrissen et~aIlllT99l . The 
comparison between the predicted and observed number 
of WR stars favors the case of an IMF slope flatter than 
2.35. 

The FUSE spectrogram of the inner part of NGC 604 
obtained through the MDRS aperture is shown is Fig- 
ure 2b (S/N~14). This smaller aperture corresponds to 
a physical size of 16x82pc 2 . The stellar line profiles are 
similar to those obtained with the LWRS aperture, but 
not exactly the same. The C III and P V line profiles 
cannot be reproduced as well as for the LWRS data, es- 
pecially in their blue wings. The models closer to the 
observed line profiles are those of 3.9-4. lMyr at 0.1 Z Q 
and 3.3-3.4Myr using 0.4 Z Q models. Interestingly, the 
MDRS spectrogram of NGC 604 corresponds better to a 
combination of a synthesized population and the spec- 
trogram of a 081 LMC star. The blue wings in P V 
and C III profiles are fitted by the single star spectro- 
gram, while the photosphcric portion cannot be fitted by 
the star, but by a modeled population. This strongly 
suggests that the number of massive stars within the 
aperture is low enough to be subject to statistical bi- 
ases on the stellar IMF, and is not well represented any- 
more by an analytical IMF. Assuming a stellar popula- 
tion of 3.3 Myr at 0.4 Z© as found previously, the con- 
tinuum slope for the MDRS spectrogram gives E(B-V)i 
0.03±0.02 if no Galactic extinction is considered. The 
flux level indicates a stellar mass of about lxl0 3 M© 
through the MDRS aperture, clearly indicating that the 
MDRS aperture does not include the whole GEHR. 

4.2. NGC595 

As NGC 604, NGC 595 con tains multipl e star 
clusters with OB stars (e.g. | Dris sen eij al.1 119931: 
Mas- Hesse fc Kunthl 119991 iMai'z-ApellanizI 120011) . The 
FUSE spectrogram of NGC 595 is presented in Figure 2c 
with S/N~13. Particularly strong P Cygni profiles are 
observed in C III and P V. As for a single evolved 
O star, the C III profile of NGG595 does not show a 
blend of photospheric+wind features as in an integrated 
population, but a single well-developed P Cygni profile. 
In fact, it appears that a synthesized stellar population 
is unable to reproduce the FUV line profiles. The FUSE 
spectrograms have then been compared to those of single 
O stars from the FUV stellar library of LavalSB and it 
reveals that an 071 LMC star is the closest match to 
the spectrogram of NGC 595 (see superimposed thick 
line spectrogram in Figure 2c). It is obvious here that 
there are not enough hot stars in NGC 595 to fit an 
analytical IMF as used in current spectral synthesis. 
Only a few stars with strong winds seem to dominate 
the line profiles. 

According to LavalSB, 071 stars appear between 
2.5 and 4.0 Myr after an instantaneous burst. At 



2.5 Myr, stars slightly brighter than 071 will proba- 
bly dominate the FUV flux. Consequently, the 071 
stars in NGC 595 would be consistent with an age 
of 3.5±0.5Myr with a metallicity close to the LMC 
(0.4 Z^). This age i s co nsistent with the works o f 
iMalumuth et~aTl (|1996[) and lMas-Hesse fc Kunthl l|1999j) . 
Assuming a standard IMF, it is still possible to roughly 
estimate parameters related to the FUV slope and flux 
level. Adopting a Galactic extinction of 0.04 (NED 4 ), 
a very low internal extinction of E(B-V)i=0.02±0.02 
is found. The stellar mass of NGC 595 is then esti- 
mated to be about lxlO 3 M with very large uncer- 
tainties. Previous works in the yisible range suggested 
a higher extinction value of 0.3 llMalum uth et al.lll996t 
iMas- Hesse fc KuntlJll1^lMaiz-Apellanizll2001l) for this 
GEHR . iMas-Hesse fc Kunthl (|1999fl and lMalumuth et alJ 
(1996) also estimated a stellar mass of 5-6xlO 3 M0, 
which is also significantly higher than the FUV result, 
but of the same order of magnitude. 

Based on LavalSB, the age and mass of NGC 595 sug- 
gest that about 10 O stars and 1 or 2 WR should be 
present in NGC 595. Howe ver, HST imaging r eveals 
larger numbers of these stars. IDrissen et, al.l 11993 1) iden- 
tified 11 WR/Of candidates and lMalumuth et all lll996ft 
estima ted the number of O stars to be ~90. IDrissen et alJ 
(1993) estimate that there are 2.5 times fewer stars be- 
tween 15 and 60 M Q in NGC 595 than NGC 604, implying 
that NGC 595 must be about 2.5 times less massive than 
NGC 604. FUV synthesis gives a factor of 5 between the 
stellar masses of the two GEHR. Recently, optical spec- 
tra from Drissen et al. (2005, in preparation) confirmed 
the presence of several WR candidates within NGC 595 
and classified them. Based on the HST/WFPC2-F170W 
archival image, the WR stars produce about 30% of the 
UV luminosity. Obviously, the observed number of WR 
stars in this object is incoherent with the FUV synthesis 
point of view. 

In an attempt to reproduce the observed FUV spec- 
trogram, simple combinations of individual hot stars are 
tested. The combinations are comprised of individual 
late O-type stars (or synthetic models) and WR stars for 
which ~ 30% of the total FUV flux comes from 1 WN6/7 
star and 4 WN7/8 stars, as classified by Drissen et al. 
(2005, in preparation). However, the resulting fits are 
poor, with the stellar combinations always giving wind 
profiles too strong in emission and having too narrow 
blue absorption. Howev er, the FUSE atlas of WR stars 
from iWillis et alJ l|2004j) revealed spectra of WR stars 
in general with spectral line profiles that are changing 
considerably from one type to another. A closer look at 
this atlas shows that HDE 269927, a WN9 type star from 
the Galaxy, display line profiles of C III and P V similar 
to stellar lines of NGC 595. Replacing the WN7/8 spec- 
tra used in the previous combinations by the spectra of 
HDE 269927 gives surprisingly good results. In fact, the 
combination of spectrograms from a 071 star (70% of 
the flux) as well as 1 WN6 and 4 WN9 stars (30% of 
the flux) reproduces well the FUSE data for NGC 595. 
This implies two things. First, it appears that the FUV 

4 The NASA Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space Administra- 
tion; http://nedwww.ipac.caltech.edu/ . 



Massive Stars in Giant H II Regions of M 33 and M 101 



5 



spectra of WR stars show line profiles that change sig- 
nificantly from one spectral type to another, and that 
probably vary with metallicity as well. Consequently, 
the few WR spectrograms currently used in the LavalSB 
spectral library are probably not very representative of 
their spectral types. Fortunately, these stars do not usu- 
ally contribute significantly to an integrated stellar pop- 
ulation and then do not really affect the synthetic spec- 
tra. Second, it seems obvious that the FUV spectra of 
NGC 595 is dominated by evolved late-type O and WN- 
late stars. However, one fundamental question remains: 
how did NGC 595 come to produce a stellar population 
enhanced in WR stars? 

The FUV synthesis of NGC 595 implies that 
F(Ha)=(1.3±0.2)xlO~ 12 ergs s" 1 c m" 2 A" 1 . Various 
values are found in the literature. IBosch et alJ (|2002f) 
obtained l.lxlO" 11 ergs s" 1 cm" 2 A" 1 , and jKennicutti 
Ijl979fl measured 8.8xl0 _12 ergs s _1 cm" 2 A" 1 . It is 
obvious that the FUV synthesis is not accurate in this 
case, and possibly also that it does not include the entire 
GEHR. 

The FUSE spectrogram of NGC 595 clearly reveals 
that a stellar population with a stellar mass of a few 
10 3 M© is too small to apply the spectral synthesis tech- 
nique, at least below 1200A. Obviously, statistical fluctu- 
ations related to a small number of massive stars are not 
well represented by an analytical IMF. A more detailed 
discussion on this subject will be given in H5.ll 

4.3. NGC 592 

Because of its fainter Ha luminosity, NGC 592 is a 
much less studied GEHR, but not less interesting. The 
observed FUV spectrogram is shown in Figure 2d, with a 
rather low S /N of 6. The FUSE aperture contai ns the en- 
tire GEHR l|Bosch et al.ll2002HKeel et al.ll2004j) . Despite 
noisy stellar lines, their profiles clearly display extended 
blue absorption wings from evolved O stars. Compar- 
ing both P V A1128.0 and C III A1175.6 lines to the 
models, it is possible to reproduce their profiles with a 
4.0±0.5Myr stellar population at Z© metallicity. Models 
at 0.4 Z© produce too weak PCygni effects in C III. The 
spectrogram is too noisy to discriminate between vari- 
ous IMF slopes. F rom Ha and H/3 narrow-band images, 
IBosch et all <|2QQ2f) estimated the age of NGC 592 to be 
more than 4.5 Myr, which is not really compatible with 
FUV line profiles displaying rela tively strong P C ygni 
features. In term of metallicity, iKeel et aTl l)2004l) in- 
terpolated a value of 0.5 Z© in [O/H], and Drissen et al. 
(2005, in preparation) estimated that 12+log[0/H]~8.4 
(i.e. 0.5 Z ) from [O III] /H/3 and [N II]/Ha line ratios. 
These values are consistent with the FUV synthesis con- 
sidering that Zq models can cover rel atively well a m etal- 
licity range from 0.4-0.5 to -1.2 Z Q l)Pellerinll2004|) . 

Using a model of 4.0 Myr at Z© and a standard IMF, 
and assuming a Galactic extinction of 0.042 (NED), an 
E(B-V)j of 0.07±0.02 is deduced from the FUV contin- 
uum slope. Once the data are corrected for extinction, 
the stellar mass deduced is (1.1±0.3) x 10 4 M Q . This 
mass is similar to that estimated for NGC 604, which 
is consistent with the fact that the stellar line profiles 
can be reproduced with a synthesis technique and an an- 
alytical IMF, contrary to NGC 595. The FUV flux level 
implies a unreddened Ha flux of (2.7±0.5)xl0 _12 ergs 



s 1 cm 2 A 1 , which is the exact value bv IBosch et alJ 
( 2002) . Other predicted parameters are reported in Ta- 
ble 3. 

4.4. NGC 588 

The FUSE spectrogram of NGC 588 is presented 
in Figure 2e, with a good S/N of 12. The FUSE 
aperture includes the entire H II region l)Bosch et alJ 
120021 IKeel et all 120041) . Models at Z© produce stel- 
lar lines definitely too deep compared to the obser- 
vations. With models at 0.4 Z© metallicity, a good 
fit can be obtained for a 3.5±0.5Myr population with 
a(IMF)<2.35. A flatter IMF tends to give better re- 
sults, but it is hard to really distinguish between var- 
ious IMF slopes because of the relatively low S/N. 
Good fits can also be obtained with 0.1 Z© models of 
4.5±1.0Myr and still with a(IMF)<2.35. In the liter- 
ature, age s of 2.8, >4.5, and 4.2 Myr are reporte d for 
NGC 588 l|Ma,s-Hesse fc Kunthl ITTiM IBosch et a.l.ll2(¥)l 
Uamet et al.l l2004 respec tively), in general a greement 
with FUV line profiles. iVilchez et all l)1988|) derived 
a precise oxygen abundance of 12+log[O/H]=8.30 (i.e. 
0.4 Z©), fa voring the models at . 4 Z^. A fla t IMF is also 
favore d bv lMas-Hesse fc Kunthl l)1999f) and Uamet et alJ 
(|2004[) obtained a(IMF)=2.37±0.16 from a star counting 
method. 

Based on the best-fit model at 0.4 Z©, a low inter- 
nal extinction of at most 0.06±0.02 is measured, which 
leads to a stellar mass of (1.3±0.6)xl0 3 M©. The mass 
is higher, (4±l)xl0 3 M©, if we consider a=2.35. De- 
pending on the extinction law used , E(B-V)j, values be- 
tween 0.11 and 0.08 are measured ijMas-Hesse fc Kunthl 
Il999t Uamet et al.ll200l . These same authors obtained 
stellar masses of 534 and 3000-5800 M©, respectively. 
The smallest value was deduced from IUE data (aper- 
ture of 10"x20"), and the largest mass is from full field 
imaging data, which explains the discrepancy. FUV 
data are in relatively good agreement with imaging data, 
which suggests that most OB stars of NGC 588 are 
within the FUSE aperture. With such a mass, the 
model predicts that F(Ha)=2.8x 10~ 12 ergs s _1 cm~ 2 
A -1 , which is in good agreement with th e value of 
2-3 xl 0~ 12 ergs s" 1 cm" 2 A - 1 measured bv iKennicuttl 
(|1979J) and IBosch et alJ l|2002ft . The best-fit model pre- 
dicts 2 W R stars in NGC 588 , which is the exact number 
found bv Uamet et alJ l)2004|) in their HST images with 
resolved stars. 

4.5. NGC588-NW 

A FUSE spectrogram has been obtained in the vicinity 
of NGC 588 (North- West). From the Digitized Sky Sur- 
vey image (see Fig. 1), this region corresponds to a rel- 
atively compact and small cluster with a faint, extended 
nebula r ring. It was first reported by iBoulesteix et all 
<ll974l their object 281 ) and also identified in the work 
of lCourtes et alJ (^87) . The ring suggests that the clus- 
ter is more evolved than those synthesized above. The 
FUSE spectrogram for this cluster is shown in Figure 2f 
(S/N~7). Diagnostic stellar lines do not display P Cygni 
profiles. Synthetic models do not reproduce well the line 
profiles. The best fit is obtained for a stellar popula- 
tion around 5-6 Myr old at 0.4 Z©, but the line profiles 
are not properly fitted. A possible alternative is a single 
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star spectrum, as was the case for NGC595. Then, a 
Galactic 09.5 III star also consistent with a population 
of 5-6 Myr, gives a better match than the model but sig- 
nificant discrepancies still exist. This age is consistent 
with the presence of the faint extended ring seen around 
NGC 588-NW in the visible range. 

To push the synthesis further, a stellar population of 
5.5 Myr at 0.4 Z Q has been considered and an extinction 
value around and a stellar mass of about lxlO 3 M 
have been roughly estimated for this cluster. This stellar 
mass is similar to the one obtained for NGC 595. The 
relatively low mass of the cluster is a logical explana- 
tion for why the synthesis technique does not work well. 
Rough estimations of predicted observable parameters 
are reported in Table 3. 

The study of NGC 588-NW gives some other clues on 
the evolution of GEHR. First, the FUSE spectrogram of 
NGC 588-NW reveals the presence of an important stel- 
lar population. However, because of its slightly greater 
age (5-6 Myr instead of -3.5 Myr for NGC 595), the neb- 
ular emission is not as strong as for NGC 595 and this 
region is consequently much less studied. It is likely that 
NGC 588-NW is representative of what NGC 595 may 
look like in —2-3 Myr. Second, the GEHR is still young 
and massive enough at this age not to have dissolved 
yet into the galaxy background. It would be interest- 
ing to search for slightly more evolved GEHR to better 
study their evolution, such as the dissipation timescale 
of clusters. This kind of cluster (i.e. still very young 
but with significantly low nebular emission) may be at 
the origin of the di ffuse UV light in starburst galaxies 
ijMeurer et all 11995)) . NGC 588-NW is consistent with 
clusters of less than 10 3 with out O-type stars, as de- 
scribed bv lChandar et all l)2005l) for the diffuse UV com- 
ponent in starbursts. A more extensive search for this 
kind of object in local galaxies could settle this issue. 

4.6. NGC 544 7 

NGC 5447 is a GEHR in the spiral galaxy M 101 
(7.4 Mpc) tha t displays several knots of star formation 
(|Bosch et al.ll2002)) . The FUSE spectrogram has a S/N 
of 12 and is shown in Figure 3a. As shown in Fig. 1, 
the FUSE aperture does not include all knots. The spec- 
trogram does not show strong wind profiles, suggesting 
that most O stars have already disappeared. Models at 
Zq metallicity produce too deep stellar lines compared 
to the observations. Models at 0.1 Z Q cannot reproduce 
both P V and C III features at the same age. The best- 
fit model is obtained at 4.5±0.5Myr with an IMF slope 
of 2.35 or flatter. This GEHR has not been extensively 
stud ied and no age has b een proposed so far for this ob- 
ject. iScowen et alJ 1)1992)) deduced an oxygen abundance 
of 8.3 in 12+log[0/H], compatible with the line depths 
of P V and C III. 

The measured FUV slope for NGC 5447 suggests that 
E(B-V),;— . From phot ographic plates and the Balmcr 
decrement . ISmitrJ $1975) estimated an extinction of 0.37, 
much larger than the FUV value. The FUV flux indicates 
a stellar mass of (1.2±0.2) x 10 5 M Q . From FUV synthe- 
sis, LavalSB predicts that F(Ha)=(5.7±0.9) x 10~ 13 ergs 
s" 1 cm" 2 A" 1 , and an EW(Ha) = 1064A f or NGC 5447. 
Using photometric data Bosch ct al. ( 2002| measured an 
Ha flux of 4.7xl0~ 12 ergs s _1 cm -2 A -1 , and lKennicuttl 



1)19791) obtained a value of 1.6x 10~ 12 ergs s _1 cuT 2 A" 1 . 
Since the GEHR is mu ch more extended t han the FUSE 
aperture (see Fig. 5 of Bos ch et all T2002 ). the factor 5- 
10 discrepancies can easily be explained. However, the 
presence of a second generation of stars contributing to 
the nebular flux but not to the FUV flux c annot be ex- 
cluded (see q5.3|) . For their knot A only, iBosch et alJ 
(2002) obtained that F(Ha)= 7.5xl0" 13 ergs s" 1 cm" 2 
A -1 , suggesting that this knot must be the princi- 
pal contributor to the FUV f lux measured with FUSE. 
iTorres-Peimbert et al.l 111989)) measured a dereddened 
equivalent width of 1096A through a 3.8" x 12.4" slit, in 
very good agreement with the FUV predictions and the 
knot A. 

4.7. NGC 5461 

NGC 5461 is a very large GEHR (>500 pc in diam- 
eter) with multiple components in M 101 iBosch et all 
2002; Kee l et alJ2004llChen. Chu fc ,Iohnsonl20051) . The 
FUSE aperture contains most of the Ha emission and 
should include most of the massive stellar content (see 
again Fig. 1). The FUSE spectrogram is shown in Fig- 
ure 3b, with a S/N of about 7. The C III feature displays 
a wind profile, implying the presence of giant and super- 
giant O-type stars. Models at Z Q do not reproduce the 
stellar line depth. The models at 0.1 Z Q give a good fit for 
a 4.0±0.2Myr stellar population and an IMF slope flat- 
ter than 2.35. A good correspondence is also obtained 
with 0.4 Z© models at 3.3±0. 2Myr, still with <a<2.35 
A mult i wavelength study from lRosa fc Benvenutil l)1994f) 
suggests an age betw een 3.0 and 4.5 Myr, compat ible 
with FUV line profiles. iLuridiana fc Peimbertl 1)200 If) de- 
duced an age between 2.5 and 3.5 Myr based on EW(H/3), 
also in general agreement with FUV line profiles. While 
the age determination me thod using EWYH/3) is not a 
recommended diagnostic l)Terlevich et alJ 12004)) . it ap- 
pears that it still gi yes good results at a such very young 
age. More recently. iChen. Chu fc Johnson! l)200 5) identi- 
fied about 12 candidate stellar clusters within NGC 5461 
of which half of them are less than 5 Myr old. The 
other clusters are probably older and do not seem to 
contribute much to the FUV flux. Abundances ranging 
from 8 .4 to 8.6 in 12+log[Q/Hl are found in the liter- 

IT 



ature llTorres-Peimbert et alJH989t IScowen et "all 119921: 



iRosa fc Benvenutil 119941 : ILuridiana fc Peimbertl 2001). 
Their observations favor the FUV synthesis models at 
0.4Z Q . 

Comparing with the modeled population of 3.3 Myr 
at 0.4 Z Q and a=1.5, the FUV continuum slope needs 
no extinction correction. The stellar mass is then 
(1.5±0.4)xl0 4 M Q . Using a standard IMF slope of 
2.35, the calc ulated stellar mass i s then (5±l)xl0 4 M Q . 
According to IRosa fc Benvenutl 1)1994)1 . the extinction 
from the Balmcr decrement is 0.23, and using an ex- 
tinction law especially designed for M101, they find a 
stellar mass of lxl0 5 M Q . According to LavalSB, the 
FUV stellar population should produce an Ha flux of 
(5±2) x 10 -13 ergs s _1 cm" 2 A" 1 , while Ha image data 
give 6.5 and 3.2xl0~ 1 2 ergs s" 1 cm~ 2 A -1 l)Bosch et alJ 
120021 iKennicutd |l979L respectively). For this popula- 
tion, the unreddened EW(Ha ) should be about 1200A. 
ITorres-Peimbert et al.l l)1989l) obtained an unreddened 
value of 1175A, in good agreement with LavalSB pre- 
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dictions. The differences between the predicted and ob- 
served extinction, nebular flux and stellar mass will be 
discussed in more details in H5.3I 

4.8. NGC5471 

NGC 5471 is another GEHR in M 101 more compact 
than NGC 5461 and NGC 5447 and may contain about 1 9 
star clusters according to lChen. Chu k Johnson) ([2005) . 
Most of the Ha emission of this H II region would 
have been included within the LWRS aperture of FUSE. 
Unfortunately, this H II region has been observed us- 
ing the MDRS aperture (4.0" x 20"), which implies that 
some OB stars are not included in the FUV spectro- 
gram presented here (see Fig. 1). Also, in the FUSE 
data, no flux has been obtained in detector 2, which af- 
fects the quality of the synthesis since the LiF2A seg- 
ment (which falls on the missing detector ) is important 
for the S/N of P V and C III lines (see iSahnow et alJ 
2000). The FUSE spectro gram is shown in Figure 3c 
with S/N=9. The C III line profile displays no obvi- 
ous wind feature. The best-fit model is obtained for a 
stellar population of 4.5±0.5Myr at 0.4Z Q . At 0.1 Z Q , 
a modeled stellar population of 3.5-4.0 Myr can also re- 
produce the observed line profiles. Because of the noise, 
a stan dard IMF as been assumed. iMas-Hesse k Kunthl 
(1999) deduced an age of 2.9 Myr for NGC 5471, which 
is too young to explain the faint P Cygni profiles ob- 
served in the FUV diagnostic line s. Oxygen abundances 
ranging from 8.0 to 8.2 (0.2-0. 3Z w |ToiTes^P^imb^r^et_alJ 
| 1989HRosa k Benvenutll994HMas-Hesse k Kunthll999t 
iBosch et al.ll2002fl are found in the literature, which is in 
good agreement with FUV synthesis. 

Adopting the 0.4 Z Q best-fit model, the comparison be- 
tween the observed and modeled continuum slopes in- 
dicate a low extinction, smaller than the uncertainties 
of 0.02. The FUV flux level sug gests a stellar mass 
of_(7±l)xlO 4 M for NGC 5471. IMas-Hesse k Kunthl 
( 1999) obtained an extinction of 0.07 in the UV range, 
which is slightly higher than the FUV extinction. The 
FUV stellar mass d educed is consistent with t he mass 
of 1.2xl0 5 M Q from lMas-Hesse k Kun th (1993), consid- 
ering the smaller aperture used with FUSE. Predictions 
reported in Table 3 are difficult to compare with the lit- 
erature because of large differences between apertures. 
However, the FUV flux prediction is always below the 
values given from la rger apertures (e.g. lKennicuttlll979T 
IBosch et al] 12002ft . iTorres-Peimbert et alJ (|1989ft mea- 
sured a dereddened EW of 575A for Ha, consistent with 
the predictions. 

4.9. NGC 5458 

NGC 5458 is an H II region smaller and fainter than 
the previous ones in M 1 01 and not much studied ex- 
cept f or its X-ray source dWang et alJl!99flt iPence et al.l 
120011 iColbert et all 12004ft . The FUSE spectrogram is 
presented in Figure 3d, and shows a S/N~10. The spec- 
trogram displays photospheric profiles without evident 
signs of winds in both P V and C III features. Sub-solar 
metallicity models produce stellar line depths too weak 
compared to the observations. The best-fit model is ob- 
tained for a 5.5-6.0 Myr old stellar population at Z Q . A 
standard IMF has been assumed since the line profiles 
are less sensitive to the IMF when evolved O stars have 



disappeared. The continuum slope indicates a low ex- 
tinction, below the uncertainties of 0.02. The flux level 
leads to a stellar mass of (l.l±O.4)xlO 5 M . Other pre- 
dicted observable parameters for NGC 5458 are reported 
in Table 3. 

5. DISCUSSION 

The massive stellar contents of several GEHR have 
been studied in detail using the FUV spectral synthesis. 
The section below focuses on the global characteristics 
of the whole sample to better understand the physics of 
GEHR in general as well as the synthesis technique in 
the FUV. 

5.1. FUV Synthesis of Small Stellar Populations 

Spectral synthesis is a powerful technique to obtain 
a good estimate of the general characteristics of young 
integrated stellar populations. However, this technique 
usually assumes that the stars follow an analytical IMF, 
and that the stars properly fill each bin of the mass func- 
tion. But how high does the mass of the population must 
be in order to be accurately described by an analytical 
IMF? The FUV is a good wavelength range to estimate 
this minimal mass for young systems. The FUV is es- 
pecially sensitive to IMF statistical fluctuations at high 
masses since only O and B stars produce many photons 
below 1200 A. Also, GEHR are very young systems and 
the disappearance of the most massive stars does not 
significantly affect the total stellar mass of the system. 

From FUV synthesis of GEHR in M 33 and M 101 (gj, 
it appears that a stellar mass greater than lxl0 3 M Q 
is needed to properly fulfill the IMF bins. As shown 
by NGC 592, NGC 604 (LWRS), and GEHR in M101, 
a stellar mass of ^lxl0 4 M© does not seem to suffer 
much of a statistical bias. However, the FUV synthesis of 
NGC 604 (MDRS), NGC 595, and NGC 588- NW reveals 
that a stellar mass closer to ~lxl0 3 M Q becomes too 
low to obtain reliable values of the age and mass of the 
star cluster because the stellar line profiles are not those 
of a standard modeled population, but those of a mix 
of a limited number of bright stars. Note that the mass 
limit needs to be higher for younger systems, where the 
dominant stars are of earlier spectral types than those 
found in a slightly older population. This is because a 
younger population needs to better fill the IMF higher 
mass bins and a more massive total stellar population is 
t hus required. 

iCerviho k Luridianal l)2004ft studied this problem from 
a theoretical point of view. The lower mass limit of a few 
10 3 M© found here for a synthesized population is fully 
consistent with their results, which suggest that the min- 
imal initial cluster mass needed for synthesis modeling 
in the U-band is about 8xl0 3 M Q for a 5 Myr popula- 
tion at 0.4 Zq. Following their calculation, this minimal 
mass can be slightly lower at shorter wavelengths like 
the FUV range. Using HST images where the stars of 
NGC 588 were resolved, Uamet et all (|2004ft obtained a 
standard IMF slope of 2.37±0.16 for NGC 588 by using 
a star counting technique and esti mated a stellar mass of 
5.8± 0.5) x 10 3 M Q , consistent with lCervino k Luridianal 
2004) and FUV synthesis. FUSE spectral synthesis of 
GEHR has clearly shown that their calculation not only 
applies to color bands, but also to stellar line profiles. 
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5.2. The Flat IMF slope of GEHR 

The stellar IMF is a matter of debate since the work 
of iSalpeterl (^555). The generally accepted slope 5 for 
the massive OB star regime at all metallicities in ev- 
ery kind of environm ent (starburst s as well as star clus- 
ters), is a=2.35 Ce.g lMassevlll998HSchaerer et"afll200a 
IGrebel fc Gallagherf l2004t IPiskunov et al.ll2004j) How- 
ever, the IMFs of GEHR derived from FUV line profiles 
seem to favor a relatively flat slope (see Table 2). Since 
FUV stellar flux is produced only by O and B stars, a 
small change in their relative numbers can affect the de- 
rived IMF slope. This result cannot be associated to a 
bias due to the FUV synthesis since several, and big- 
ger, young populations have been studied with t he same 
techn ique and did not show such a flat slope l)Pellerinl 
l20?)l . 

Some hypotheses could physically explain a flat IMF 
in the FUV range. One hypothesis is that B-type stars 
could still be more extinguished by dust than earlier type 
stars. If so, it would then be more difficult to see them 
in the FUV, producing an artificially flatter IMF. How- 
ever, the ex tinction values of individ ual stars in NGC 604 
obtained bv lBruhweiler et "aTl <)2003j) do not show a signif- 
icant correlation with the spectral type, suggesting that 
B stars are not systematically more extinguished than O 
stars. 

Another more plausible possibility is that the massive 
stars fill the IMF high mass bins relatively well, but not 
perfectly. If some spectral types have slightly deviant 
numbers from the analytical IMF, it will slightly change 
the integrated stellar line profiles in the same direction as 
NGC 595 or NGC 604-MDRS, i.e. by accentuating the in- 
tegrated wind profiles. Since a flatter IMF also produces 
more pronounced P Cygni profiles, it would be hard to 
differentiate the two cases. Consequently, even if the 
population synthesis gives reliable and precise results on 
most physical parameters of the population (age, mass, 
metallicity, colors, fluxes) for a >lxlO 3 M0 population 
ft |5.1(l . it appears that the stellar IMF slope derived from 
the FUV line profiles is a sign of a non-perfect filling of 
the IMF high mass bins. This last possibility is sup- 
ported b y the IMF obtained from the star counting tech- 
nique of l.lamet et aTl IpOfil on NGC 588. They derived 
a standard IMF slope, but their IMF histogram clearly 
shows that some mass bins, especially at higher masses, 
are clearly deviant from the analytical slope. 

5.3. A second generation of stars in NGC 5461 

The spectral synthesis of FUSE data on NGC 5461 has 
predicted much lower values for the Ha flux (factor of 
10), the stellar mass (factor of 2 to 10), and the extinction 
than has been reported in the literature. These discrep- 
ancies are hard to explain since most of the Ha emission 
is included within the FUSE aperture. One plausible ex- 
planation is the presence of a second generation of stars 
in N GC 5461, like the on e observed in the LMC Cluster 
Nil (Wal born fc Parkerlll99l . In the case of the star- 
forming region Nil, the central region is composed of a 



3.5 Myr stellar population which dominates the UV flux. 
A surrounding nebulae is excited by a younger genera- 
tion of stars which is not obse rved at short wavelength s 
because it is heavily reddened l|Walborn &: Parkerl l992). 

The presence of a second generation of stars in 
NGC 5461, younger and consequently more extinguished 
than the first one, could explain the larger extinction 
deduced at longer wavelengths, the stellar mass discrep- 
ancy as well as the excess in nebular emission. The sec- 
ond generation cluster must then be relatively massive to 
explain the large differences in flux and mass. It is not 
excluded that younger stars from different clusters are 
present rather than a single second generation. Although 
there is no proof of such a population within NGC 5461, 
this H II region is a good candidate to host very massive 
stars, younger than those actually detected with FUSE. 

It is also possible that younger stars are present within 
other GEHR studied here. Unfortunately, because the 
FUSE aperture does not always include the whole sys- 
tem, it is impossible to confirm here if the difference 
between the predicted and observed Ha fluxes comes 
from a second generation or not, as it is the case for 
NGC 604 for example . Con sidering the detailed work of 
Mafz-ADcllaniz et al. (2004) on the attenuation maps of 
NGC 604, the differences in Ha fluxes and stellar masses 
in GEHR, including NGC 5461, might also be due, at 
least partly, to the complexity of the gas and dust spa- 
tial distribution. 

6. SUMMARY 

The evolutionary spectral synthesis technique in the 
FUV has been used to study the massive stellar content 
of nine GEHR in M33 and M 101. Stellar masses, inter- 
nal extinctions, and ages have been obtained for most of 
them. The comparison of the FUV synthesis results with 
values obtained from previous available works in various 
wavelength ranges has shown that the technique is reli- 
able in most cases. The comparison of the GEHR with 
each other has confirmed observationally that the syn- 
thesis technique must be applied to stellar populations 
of at least a few 10 3 M Q in the FUV to avoid statisti- 
cal fluctuations of the high mass end of the stellar IMF. 
It has also revealed that a flat IMF slope is apparently 
favored for GEHR in the FUV, which is likely the first 
apparent effect of statistical fluctuations of the IMF for 
low mass populations. FUV data suggests that giant H II 
regions reach their maximum nebular luminosity around 
3.0-3.5 Myr, coincident with the WR phase. Finally, the 
H II region NGC 5461 in M101 is a good candidate to host 
a second generation of stars more extinguished than, and 
formed after the cluster actually detected with FUSE. 
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5 A slope of 2.35 is traditionally called a Salpeter slope. However, 
this terminology is not appropriate for stellar masses covered by 
FUSE since the work of Salpeter] 119551) applies to a lower mass 
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Fig. 1. — FUSE aperture locations. All images are from lBosch et alj i200l) in Ha (narrow-band filter) except for NGC 588 and NGC 5458 
from DSS. Squares represent the LWRS aperture (30"x30") and the rectangles represent the MDRS aperture (4"x20"). North is up, east 
is left. 
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Fig. 2. — FUV spectrograms of five GEHR in M33 superimposed on the LavalSB best-fit model of a stellar population, a) NGC 604 
through the LWRS aperture and a 3.3 Myr model at 0.4 Z with a(IMF)=1.5. b) NGC 604 through the MDRS aperture and the same 
model as previously, c) NGC 595 with the spectrogram of an 071 LMC star, d) NGC 592 and a 4.0 Myr model at Z with a(IMF)=2.35. 
e) NGC 588 and a 3.5 Myr model at 0.4 Z Q with q(IMF)=1.5. f) NGC588-NW and the spectrogram of an 09.5111 Galactic star. Thin 
lines: FUSE spectra. Thick lines: synthetic models from LavalSB or stellar spectrograms. See text and Table 2 for details on the best-fit 
models. 
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Fig. 3. — FUV spectrograms of four GEHR in M 101 superimposed on the LavalSB best-fit model of a stellar population, a) NGC 5447 
and a 4.5 Myr model at 0.4 Z with a=2.35. b) NGC 5461 and a 3.3 Myr model at 0.4 Z Q with o=1.5. c) NGC 5471 and a 4.5 Myr model 
at 0.4 Z with a=2.35. d) NGC 5458 and a 5.5 Myr model at Z Q with a=2.35. Thin lines: FUSE spectra. Thick lines: synthetic models 
from LavalSB. See text and Table 2 for details on the best-fit models. 



TABLE 1 

FUSE Observational Parameters of GEHR in M33 and M101 



Hn Region 


a(J2000) 

h m s 


<5(J2000) 
/ // 


FUSE 


Aperture 3, 


Vhel 


F(1150) b 






Program 




[km s- 1 ] 


[ergs s _1 cm~ 2 A -1 


] 


NGC 588 


01 32 45.5 


+30 38 55 


A08604 


LWRS 


-174 


(10±l)xl0- 


14 


NGC 588-NW 


01 32 37.7 


+30 40 06 


A06105 


LWRS 




(5.1±0.5)xl0" 


14 


NGC 592 


01 33 12.3 


+30 38 49 


A08602 


LWRS 


-162 


(10±l)xl0- 


14 


NGC 595 


01 33 33.6 


+30 41 32 


A08603 


LWRS 


-174 


(9.2±0.9)xl0" 


14 




01 33 33.6 


+30 41 32 


B01801 


LWRS 




(10±l)xl0- 


14 


NGC 604 


01 34 32.5 


+30 47 04 


A08601 


LWRS 


-226 


(100±10)xl0- 


14 




01 34 32.4 


+30 47 04 


B01802 


MDRS 




(16+2) xl0~ 


14 


NGC 5447 


14 02 28.6 


+54 16 11 


B01803 


LWRS 


152 


(8.4+0.8) x 10" 


14 


NGC 5458 


14 03 12.6 


+54 17 55 


C07001 


LWRS 


241 


(3.8±0.5)xl0" 


14 


NGC 5461 


14 03 41.3 


+54 19 05 


A08605 


LWRS 


298 


(6.3+0. 7)xl0" 


14 


NGC 5471 


14 04 28.7 


+54 23 49 


B01805 


MDRS 


297 


(4.2+0. 6)xl0- 14 





a The LWRS aperture is 30"x30" and MDRS aperture is 4"x20" 
k Not corrected for extinction. 
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TABLE 2 

FUV Synthesis Parameters of GEHR in M33 and M 101 



Hll Region 


Age 




a (IMF) 


E(B-V)< 


Stellar Mass 




[Myr] 






±0.02 


[M ] 


Best-Fit Parameters 










NGC 588 


3.5±0.5 


0.4 


<2.35 


0.06 


(1.3±0.6)xl0 3 


NGC 588- NW 


5-6 


0.4 


- 


~0 


~lxl0 3 


NGC 592 


4.0±0.5 


1 


2.35 


0.07 


(1.0±0.3)xl0 4 


NGC 595 


3.5±0.5 


0.4 


- 


0.02 


~lxl0 3 


NGC 604-LWRS 


3.3±0.1 


0.4 


1.50 


0.03 


(7±2)xl0 3 


NGC 604-MDRS 


3.3±0.2 


0.4 


1.50 


0.03 


(1.0±0.3)xl0 3 


NGC 5447 


4.5±0.5 


0.4 


<2.35 


~0 


(1.2±0.2)xl0 5 


NGC 5458 


5.5±0.5 


1 


2.35 


~0 


(l.l±0.4)xl0 5 


NGC 5461 


3.3±0.2 


0.4 


1.50 


~0 


(1.5±0.4)xl0 4 


NGC 5471 


4.5±0.5 


0.4 


2.35 


~0 


(7±l)xl0 4 


Other Good-Fit Parameters 








NGC 588 


4.5±1.0 


0.1 


<2.35 






NGC 604-LWRS 


3.9±0.1 


0.1 


1.50 


0.03 


(l.l±0.2)xl0 4 


NGC 604-LWRS 


3.5±0.3 


0.1 


2.35 


0.03 


(2.5±0.3)xl0 4 


NGC 604-LWRS 


3.0±0.3 


0.4 


2.35 


0.03 


(2.2±0.5)xl0 4 


NGC 5461 


4.0±0.5 


0.1 


1.50 


~0 


(2.2±0.6)xl0 4 


NGC 5471 


3.5-4.0 


0.1 


2.35 


~0 



TABLE 3 

Predicted Parameters for GEHR from FUV Synthesis 



Hii Region 


#0 #WR F(H«) a F(5500) a 


F(Ha) a 


Ref 




Predicted 


Observed 





NGC 588 


15 


2 


2.8x10" 


-12 


4x10" 


-15 


2-3x10" 


-12 


1, 2 


NGC 588-NW 


4 





2x10" 


-13 


1x10" 


-IS 








NGC 592 


40 


1 


2.7x10" 


-12 


1x10" 


-14 


2.7x10" 


-12 


2 


NGC 595 


10 


1 


1.3x10" 


-12 


2x10" 


-15 


8.8-11x10" 


-12 


2, 1 


NGC 604-LWRS 


90 


9 


2x10" 


-11 


3x10" 


-12 


3.3-4.0x10" 


-11 


3, 2 


NGC 604-MDRS 


12 


2 


2x10" 


-12 


3x10" 


-13 








NGC 5447 


460 


34 


5.7x10" 


-13 


9.2x10" 


-IS 


1.6-4.7x10" 


-12 


1) 2 


NGC 5458 


150 


46 


1.7x10" 


-13 


6x10" 


-16 








NGC 5461 


175 


18 


5x10" 


-13 


7x10" 


-16 


3.2-6.5x10" 


-12 


1, 2 


NGC 5471 


290 


20 


3.5x10" 


-13 


6x10" 


-16 


3.5-6.2x10" 


-12 


1) 2 
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